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Dynamic  ductile  fracture  of  a  two-phase  material  has  been  investigated 
both  experimentally  and  analytically.  This  study  consists  of  two  stages.  The 
first  stage  is  "a  model  study"  where  the  idealized  two-phase  material,  single 
crystal  Cu-SiC>2  use<j  for  both  experimental  and  analytical  studies.  For  the 
analytical  study,  the  growth  of  voids  was  focused  on,  and  two  analytical  models 
(static  and  dynamic  void  growth  models)  were  constructed  to  simulate  the 
simultaneous  growth  of  voids  which  were  nucleated  from  secondary  particles. 
The  experimental  work  was  conducted  by  using  plate  impact  test,  which  generated 

the  distribution  of  voids  in  a  Cu-Si02  plate,  and  the  results  of  the  void 
density  f  were  compared  with  those  predicted  by  the  static  and  dynamic  growth 
models.  For  smaller  strain  (or  strain-rates),  the  static  void  growth  model  can 
predict  the  experimental  results  reasonably  well,  but  it  fails  to  predict  for 
the  larger  strain- rates.  While  the  dynamic  void  growth  model  can  predict  the 
experimental  results  reasonably  well  for  the  entire  range  of  3train-rates. 
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TEM  study  conducted  on  the  as-impacted  Cu-Si02  specimens  has  revealed  that 
dislocation  cells  were  formed  at  higher  strain-rates  with  the  average  cell 
diameter  being  0.5  urn,  and  also  that  many  microcrystallites  ~io-50  nm  in 
-  diameter  were  formed  adjacent  to  the  growing  voids. 
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The  second  3tage  of  this  study  is  the  experimental  work  on  more  realistic 
two-phase  material  systems,  Al-Si  alloy  and  SiC  whisker/2124  A1  composites 
which  were  subjected  to  uniaxial  tension  by  3plit  Hopkinson  bar  (SHB)  with 
strain-rates  10~3 -lo^sec.  Then,  the  void  density  f  was  measured  along  the 
axis  of  the  fractured  specimen.  The  measured  f  was  compared  with  the 
analytical  results,  resulting  in  a  reasonably  good  agreement  between  them. 
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TEM  Photo  if  .single  crystal  Cu--3iO,  pla-e  before  impact  test. 

Schematic  view  of  plate  impact  test  apparatus. 

SEM  picture  showing  voids  produced  across  the  thickness  at 
various  impact  velocities:  (a)  V  =  no  m/3;  <  b>  V„  =  160  m/'s; 

<c>  V,  =  208  m/s. 

Distribution  of  the  volume  fractions  of  void3,  f,  across  the 
plate  thickness  at  various  impact  velocities  The  error  of 

the  peak  value  is  indicated  by  arrow. 

Laarangian  diaaram  at  impact  velocity  160  m/s. 

Void  nucleated  at  a  group  of  SiO,  particles  after  strainina  at  e, 
=  8,10-,-.  “  A 

TEM  photos  of  samples  tested  at  e,  =  la)  5.5  x  10*  s'1.  lb)  8  x 

A 

10*  s  1  and  ici  10.4  x  10s*  s'1  showing  thickening  of  the  cell 
walls  with  increasing  strain  rate. 

Sample  strained  at  *  i.04  x  io5  s"1  showino  microcrystallites 

adiacent  to  a  large  void. 

Split  Hopkinson  bar  apparatus  <a>  and  its  fixture  for  tensile 
specimen  'b1. 

Stress- strain  curves  of  Al-Si  alloy  at  various  strain-rates. 
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Stress-strain  curves  of  3iCw/2124  A1  at  ‘e  =  60o/sec  with  various 
volume  fraction  of  SiC  whisker. 

Figure  13  Void  density  measured  along  the  axis  of  the  tensile  SiCw/2i24  A1 
specimens  with  various  volume  fractions  of  SiC  whisker. 

Figure  14  Analytical  model  for  the  growth  of  voids. 

Figure  15  The  volume  fraction  of  voids  calculated  by  the  static  void  growth 

model  (f)  vs  the  applied  strain  (4ftt)  for  the  uniaxial  strain 
mode. 

Figure  16  The  volume  fraction  of  voids  (f)  calculated  by  the  dynamic  void 
growth  model.  The  results  for  fQ  =  0<005  and  0.0008  are  plotted 
by  solid  and  dashed  curves,  respectively  and  the  experimental 
results  by  open  symbols  with  the  error  arrow. 

Table  1  Experimental  results  of  the  plate  impact  test  on  single  crystal 

Cu-SiOz. 

Table  2  Comparison  between  the  analytical  results  based  on  the  static 

void  growth  model  and  experimental  results. 
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II.  Statement  of  the  Problem 

Most  commercial  metal  alloys  contain  secondary  particles  (or  fillers).  If 
the  matrix  is  ductile,  the  failure  mode  of  the  metal  is  ductile  fracture,  and 
its  mechanism  under  quasi-static  deformation  has  been  well  studied  Cl-133. 
Based  on  these  studies  it  was  found  that  voids  are  nucleated  at  secondary 
particles  and  they  grow  continuously  with  the  applied  strain  (stress),  and 
ultimately  coalesce,  resulting  in  final  fracture  of  the  metal. 

Unlike  the  case  of  quasi- static  ductile  fracture,  only  a  limited  amount  of 
work  has  been  carried  out  on  the  micro-mechanical  aspects  of  dynamic  ductile 
fracture  of  metals  C14-173.  The  previous  experimental  works  on  the  dynamic 
ductile  fracture  of  metals  are  focused  on  several  metal  systems,  high  purity 
aluminum  and  copper  C14,153  and  metal  alloys  Q4,163  which  are  plate  impacted. 
As  to  the  analytical  work  on  the  dynamic  ductile  fracture,  the  mechanism  of  the 
initiation  and  growth  of  voids  has  been  focused  on  C14,163.  The  essence  of  the 
above  analysis  is  based  on  the  assumption  that  the  matrix  metal  is  linear 
viscous  fluid  C183.  On  the  other  hand,  Glennie  C173  treated  the  void  growth  in 
a  perfectly  rigid  plastic  material  by  using  the  Hill's  variational  method  C193 
and  modifying  the  static  void  growth  model  by  Rice  and  Tracey  C203. 

The  above  brief  review  on  the  previous  works  on  the  dynamic  ductile 
fracture  of  metals  had  led  us  to  believe  that  the  following  subjects  need  to  be 
focused  on: 


the  mechanism  of  dynamic  ductile  fracture  of  an  idealized  two-phase 
material  system  (single  crystal  Cu-SiO^. 

the  modeling  of  the  dynamic  void  growth  in  a  non-linear  viscous 
metal. 

the  dislocation  morphology  in  an  as-impacted  two-phase  material, 
the  mechanism  of  dynamic  ductile  fracture  of  a  new  class  of  two-phase 
material  including  metal  matrix  composites  lAl-Si  alloy,  SiCp  and 


s-i-cw/Al  composites). 
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In  order  to  study  the  above  subjects,  we  have  conducted  the  following  tasks: 
Experimental  Work 


ft 


% 


v. 


(i)  Plate  impact  test  on  single  crystal  Cu-SiO^  wit^  several 
strain-rates. 

(ii)  Split  Hopkinson  bar  (SHB)  test  on  Al-Si  alloy  and  SiCw/^ 
composite. 

(iii) Drop  weight  impact  test  on  SiCp/fti  composite. 

(iv)  Metallographic  study  on  the  as-impacted  specimen. 

Analytical  Work 

(i)  Construction  of  the  static  void  growth  model. 

(ii)  Construction  of  the  dynamic  void  growth  model. 

III.  Summary  of  the  Most  Important  Findings 

Both  experimental  and  analytical  works  have  led  to  a  number  of  important 
,  findings,  which  will  be  given  below. 
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vl.  Experimental  Findings 


Single  crystal  Cu-SiO^  (Fig.  1)  plates  were  plate  impacted  at  various 
striking  velocities,  Vg  =  no,  160  and  208  m/sec  (see  Fig. 2  for  the  plate 
impact  test  apparatus).  The  impacted  specimens  were  sectioned  and 
examined  by  SEM  (Fig. 3),  and  the  void  density  f  across  the  plate  thickness 
was  measured  (Fig.  4).  The  experimental  results  are  summarized  in  Table 
1.  The  as-impacted  specimen  was  also  examined  by  TEM.  The  SEM  and  TEM 
studies  have  revealed  the  following  important  findings  C22D: 


v 
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a.  Void  density  f  increases  with  the  duration  of  the  tensile  stress  wave 
'At,  see  tig.  5  for  its  definition)  and  Vg>  This  observation  is 
consistent  with  that  by  Seaman  et  al  C14D. 


b.  Only  a  fraction  of  secondary  particle  Si02  have  nucleated  voids  ( see 
Fig.  6),  which  is  in  marked  contrast  with  the  case  of  the 
quasi-static  loading  where  virtually  all  the  secondary  particles 

<Si02)  nucleate  voids  C21J. 

c.  TEM  study  revealed  the  formation  of  dislocation  cell  structures  in 
the  as-impacted  single  crystal  Cu-Si02  wj.th  the  average  cell  diameter 
being  constant  (0.5  my)  for  the  range  of  strain-rates  investigated 

(104  -  105/ sec;  (see  Fig.  7),  and  also  that  microcrystallites  were 

formed  near  a  large  void  (Fig.  8).  The  formation  of 

microcrystallites  is  attributed  to  either  dynamic  recrystallization 
or  recrystallization  due  to  adiabatic  heating. 

In  addition  to  the  idealized  two-phase  material,  we  have  also  studied 
the  dynamic  ductile  fracture  of  ordinary  two-phase  metals.  Al-Si  alloy 
and  SiC  whisker/2124  aluminum  (SiC^^)  composite  by  subjecting  them  to 
uniaxial  tension  (Split  Hopkinson  Bar  (SHB),  3ee  Fig.  9)  at  various  high 
strain-rates  up  to  1500/sec.  The  dynamic  stress-strain  curves  of  these 
two-phase  materials  were  obtained.  Then  the  fractured  specimens  were 
sectioned  along  the  specimen  axis  for  the  examination  of  the  initiation 
and  growth  of  voids  by  SEM.  The  above  experimental  study  on  the  ordinary 
two-phase  materials  has  led  to  the  following  findings: 

d.  Both  Al-Si  alloy  and  SiC^^  composite  exhibited  strong  strain-rate 
sensitivity  (see  Fig.  10  for  Al-Si  alloy  and  Fig.  11  for  SiCw/Ai  wdth 
volume  fraction  of  SiC  whisker  beinq  15%). 

e.  The  higher  the  volume  fraction  of  secondary  filler  (f  >r  the  larger 
the  flow  stress  (Fig.  12)  and  also  the  higher  the  void  density  f 
(Fig.  13)  become. 

SiC  particulate/6061  A1  (SiCp/Al)  composite  was  impacted  by  drop 
weiqht  test  and  then  subjected  to  three-point  bendinq  test  to  obtain  the 
residual  strength  1293.  This  drop  weight  test  generates  lower 
3train-rate,  10/sec.  This  study  has  yielded  the  following  findings: 

f.  Voids  were  not  observed  in  the  as-impacted  specimen. 


The  residuai  fracture  energy  or  the  as -impacted  composite  decreases 
with  the  increase  in  the  impact  eneruy  and  the  rate  of  the  reduction 
in  the  fracture  energy  becomes  larger  for  smaller  volume  fraction  of 
particles. 

Table  1.  Experimental  Results  L22J 


Symbols 

Case  1 

Case  2 

Case  3 

Striking  velocity  of  projectile  (ra/s) 

V 

s 

no. 

160. 

208. 

Hugoniot  elastic  limit  (GPa) 

ae 

0.164 

0.164 

0.164 

Elastic  particle  velocity  (mm/us) 

u 

e 

0.00415 

0.00415 

0.00415 

Elastic  stress  wave  velocity  (mm/ps) 

U 

e 

4.46 

4.46 

4.46 

Plastic  particle  velocity  (mm/us) 

u 

P 

0.055 

0.080 

0.104 

Plastic  stress  wave  velocity  (mm/ps) 

U 

p 

4.0010 

4.0393 

4.0747 

Initial  density  of  Cu  (mm/ys) 

po 

8.93 

8.93 

8.93 

Density  of  Cu  behind  the  plastic 

shock  front 

pi 

9.01 

9.06 

9.10 

Duration  (ys) 

At 

0.186 

0.190 

0.195 

Strain-rate  (104/s) 

5.5 

8.0 

10.4 

Applied  strain  (1CT2) 

' 

eAat 

1.023 

1.52 

2.028 

Analytical  Findings 

In  order  to  predict  the  void  density  observed  in  the  as-plate 

impacted  Cu-Si02  and  also  the  as-SHB  tested  ordinary  two-phase  materials, 
we  have  constructed  two  analytical  models,  static  void  growth  model 
Cll.13,243  and  dynamic  void  growth  model  C22.26D.  The  semi-detailed 
statement  of  the  static  and  dynamic  void  growth  model  is  given  in 
Appendices.  It  was  assumed  in  our  model  that  voids  of  the  same  size  will 
grow  simultaneously  under  applied  strain,  where  eA  and  t  are  applied 

strain-rate  and  time,  respectively  (Fig.  14).  In  our  models,  we  focus  on 
a  representative  void  < the  solid  circle  in  Fig.  14)  which  is  surrounded  bv 


a  mixture  of  voids  ind  incompressible  matrix  metal  and  have  calculated  the 
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growth  of  the  representative  void.  Thus,  the  surrounding  media  becomes 
compressible  a3  far  as  the  volume  fraction  of  voids  (f)  >  0.  In  this 

respect,  the  interaction  between  growing  voids  can  be  taken  into  account, 
whereas  the  other  models  cannot  account  for  this  interaction  since  they 
focus  on  a  single  void  surrounded  by  incompressible  matrix.  The 
analytical  study  has  led  to  the  following  findings. 

a.  A  comparison  between  the  experimental  results  of  the  peak  values  of  f 
(Fig.  4)  and  the  analytical  results  by  our  static  void  growth  model 
C243  revealed  that  the  static  void  growth  model  can  predict  the  void 
density  f  reasonably  well  for  smaller  but  it  fails  for  larger 

e^t  (Fig.  15).  f0  in  Fig.  15  denotes  the  initial  volume  fraction  of 
voids  and  should  be  less  than  or  equal  to  0.005  which  is  the  volume 

fraction  of  the  secondary  particle  Si02.  fQ  wa3  found  to  be  about 
0.0008  C223. 

b.  Our  dynamic  void  growth  model  11261  can  account  for  the  interaction 
between  growing  voids,  material  non-linearity  (1  <  n  <  »)  and  inertia 
effect.  Besides,  this  model  can  yield  the  existing  solutions  as 
special  cases,  i.e.,  the  solutions  of  the  dynamic  void  growth  by 
Poritsky  for  linear  viscous  metal  C183,  and  by  Gennie  for  perfectly 
rigid  plastic  material  (n  =  <»)  C173,  and  those  of  the  static  void 
growth  by  Budiansky  et  al  C103.  The  analytical  results  by  our 
dynamic  void  growth  model  can  predict  the  experimental  results  of  the 
peak  value  of  f  in  the  as-plate  impacted  Cu-Si02  well  (see  Fig.  16). 
It  can  be  concluded  from  Fig.  16  that  the  matrix  metal  (Cu)  exhibits 
a  weak  non-linearity,  n=1.5  (see  the  definition  of  n  in  Appendix  B). 

c.  Another  important  finding  in  the  dynamic  void  growth  by  our  model 
C263  is  that  the  non-dimensional  parameter  y  defined  by  strain-rate, 
initial  radius  of  the  void  and  material  parameters  (see  the 

definition  of  yQ  Appendix  B)  play  a  dominant  role  on  the  dynamic 
void  growth.  Namely,  when  yQ  -j_s  iarq9,  the  inertia  effect  is  to 
reduce  the  growth  rate  at  the  initial  stage  of  void  growth.  If  yQ  i3 
small,  the  dynamic  void  qrowth  can  be  approximated  by  the  static 
growth  model.  The  material  non-linearity  has  a  strong  effect  on  the 


7 


void  growth  for  smaller  yQ  ( quasi  static  growth),  but  a  weak  effect 
on  the  void  growth  for  larger  yQ  (dynamic  growth.. 

The  void  density  measured  in  the  as-SHB  tested  ordinary  two-phase 
material  was  compared  with  that  predicted  by  the  static  void  growth 
model,  resulting  in  a  reasonably  good  agreement  though  the  analytical 
prediciton  fails  to  explain  the  trend  of  the  increase  in  f  as  e 
increases  (see  Table  2)  C303.  Then  the  dynamic  void  growth  model  was 
applied  to  this  case,  leading  to  the  better  agreement  between  the 
experimental  and  analytical  results.  In  using  the  dynamic  model,  the 
ratio  of  the  mean  stress  to  the  effective  stress  was  properly 
estimated  by  use  using  the  results  by  Needleman  C273. 


Table  2  Comparison  between  the  experimental  and 
analytical  results  C303 


strain-rate,  e 


Location  from  the 
fracture  section 


Void  density 
measured 

_ _ _ i _ 

2  mm 

0  mm 

0.0013 

0.005 

0.009 

0.005 

0.025 

0.035 

Void  density 
predicted  by  the  static 
void  growth  model 


2 
2 

0.0021 


0.012 

0.012 

0.013 


free  surface 


mid-plone 
Thickness  direction 


impact  surface 


o.2Pm 
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VII.  Appendices 


>< 


Appendix  A: 


Static  Void  Growth  Model 


It  is  assumed  in  our  theoretical  model  that  the  effect  of  inertia  can  be 
neglected  at  the  early  stage  of  void  growth,  the  voids  grow  simultaneously  and 
the  matrix  is  incompressible  viscous  material  at  high  strain-rate:  a  =  2uo‘e 
where  a  is  the  stress,  e  is  the  strain-rate  and  is  the  viscosity.  Here  the 
two-phase  material  consists  of  a  ductile  matrix  obeying  a  =  2uQe  and  spherical 
particles  which  nucleate  spherical  voids.  After  voids  are  nucleated,  they  grow 
under  the  applied  strain-rate  which  is  assumed  to  be  constant.  Then  the 
composite  that  consists  of  the  incompressible  matrix  and  voids  becomes 
compressible  due  to  the  existence  of  the  voids.  In  order  to  solve  the  problem 
of  void  growth,  we  focus  on  a  representative  void  that  is  surrounded  by  the 

composite  media  subjected  to  eft  along  the  x3-axis  as  shown  in  Fig.  14.  The 
coefficient  E,  in  Fig.  14  denotes  either  the  case  of  uniaxial  tension  (£  =  -v), 
or  uniaxial  strain  =  0)  or  that  of  triaxial  stress  (£;  =  1).  The  composite 
has  two  constants,  the  effective  viscosity  (pi)  and  the  effective  expansion 
coefficient  (k),  and  they  are  given  by  C11J. 


M  -  UQ  U  +  3(1_f) 


K  = 


(1-f) 


(1) 


(2) 


is  the  viscosity  of  the  incompressible  matrix  material  and  f  is  the 


where 

volume  fraction  of  voids.  For  later  convenience,  the  domain  of  the 
representative  void  is  denoted  by  fl  and  that  of  the  whole  composite  by  D,  then 
that  of  the  surrounding  media  becomes  D-Q. 
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Applying  Eshelby's  equivalent  inclusion  method  C251,  we  have  in  SI 


a  6, 


.0 


ij(ekk  +  ekk  -  ekk'  +  <^<e°ij  +  -  ’e^)  =  0 
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(3) 
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where  \  =  k  -  2p/3,  u  and  k  are  defined  by  eqs.  (1)  and  (2),  x.  .  ' 

n  0 i]  is  the 

Kronecker's  delta,  the  super  dots  denote  the  material  time  derivative,  e.°  . 

•*  1S 
the  applied  strain-rate,  e^  the  eigenstrain-rate  that  has  some  value  in  Si, 

but  vanishes  in  D  -  SI,  and  e^  the  disturbance  of  the  strain-rate  due  to  the 

existence  of  £1  and  is  related  to  6*^  a3  [251 

eij  =  sijkJiekil  <4) 

where  j,3  the  EShelby's  tensor  and  its  explicit  expressions  for  SI  being 

sphere  and  prolate  ellipsoid  are  given  in  Appendix  C.  By  setting  i  =  j  in  eq. 
(3)  and  noting  that  3x  +  2\x  =  K  t  0  due  to  voids,  we  obtain 

ekk  +  ekk  ~  ekk  =  0  (5) 

where  the  repeated  index  is  to  be  summed  over  1,  2  and  3.  Also  by  setting  i  = 
j  =  3  in  eq.  (3),  we  have 

^33ilell  +  (S3333  ■  1)e33  =  “®A  (6) 

where  we  have  used  the  fact  that  the  system  of  Fig.  14  gives  rise  to  transverse 
isotropy  and  hence  only  independent  components  are  those  with  ij  =  11  (or  22) 
and  33.  The  shape  of  ft  will  become  sphere  for  £  =  1  and  prolate  ellipsoid  for 
E,  =  -v  and  0  under  the  applied  strain-rate  Let  U3  denote  the  axes  of  the 

major  and  minor  of  the  prolate  ellipsoid  by  c  and  a.  Then  the  growth  rate  of 
the  axes  are  given  by 

a  .o 

a  =  en  =  eii  (7) 

c  .o 

c  =  e33  +  e33  (8) 

After  having  solved  tor  e^  and  e'3  *  tn  eq3.  <5>  an(j  (6)r  then  using  eqs.  (4), 
(7)  and  (8),  we  arrive  at  for  triaxial  stress  (£  =  1), 


for  uniaxial  strain  (£  =  0), 


4  =  c  S1133(S1111  *  S1122)  *  (1  ~  Sllll  ~  S1122)S1133  y& 
a  (1  "  Sllll  ~  s1122)(1  "  s3333>  ‘  2S113333311  A 

c  =  Q  +  ZS1133S3311  1  (1  ~  Sllll  S1122)S3333 _  u 

c  (1  “  Sllll  ‘  s1122 ’ (1  "  S3333>  '  2S1133S3311  A 

and  for  uniaxial  stress  (£  =  -v), 

a.,.  v+  CS1133~v(1~S3333)3(SlllltS1122)'t'(1~Sllll"S1122~2vS3311)S1133K 
a  V+  C(1'S1111_S1122)(1“S3333)-2S1133S3311} 


^_c1t2S3311<:S1133~v(  1~S3333 }  }+a~Smi'S1122"2vS3311}S3333_-|£ 
c  C(1‘Sllll-Sii22><l-s3333)_2S1133s331l5 

Once  a  and  c  are  computed  at  current  time,  then  the  current  volume  fraction  of 
void,  f,  can  be  obtained  from 


(14) 


In  eqs.  (12)  and  (13),  v  is  defined  by 


Dynamic  Void  Growth  Model 


Here,  we  use  the  same  theoretical  model  as  that  for  the  static  void  growth 
(Fig.  14)  except  that  the  non-linearity  of  the  material  is  introduced  in  the 
present  model. 

The  infinite  body  D  is  subjected  to  the  applied  strain-rate  at  far  field, 

•  00 

ell  =  £22  =  and  e33  =  ®A  as  shown  in  Fig.  14.  The  strain-rate  field  is 
axi- symmetric  with  respect  to  the  x-j-axis  an(j  the  corresponding  stresses  are 
also  axi- symmetric  with  respect  to  x^_axts  an(j  can  jje  obtained  explicitly  in 
terms  of  applied  strain-rate  if  the  constitutive  equation  of  the  material  is 
known. 

The  non-linear  material  considered  is  isotropic  and  compressible  with 
power  law  dependence  of  the  strain-rate  on  stress.  Under  simple  tension,  the 
total  strain-rate  <e)  is  related  to  the  stress  (a)  by 


e  =  e 


(— ,n 
o 'a  ' 
o 


(1) 


-  Mhere  %  and  £q  are  the 

reference  stress  and 

strain-rate,  respectively. 

The 

vrate  hardening  exponent  n 

ranges  from  1  to  <». 

For 

the  multiaxial 

case. 

the 

corresponding  deviatoric 

stress  (sij}  can 

be 

related  to  the 

deviatoric 

strain-rate  (e*^)  a3 

ij 


(2) 


where  <*e  effective  stress  defined  by  ae  =  ( |3ijSij  )1/2  [28D,  and  the 
repeated  indices  are  to  be  summed  (this  rule  will  be  used  throughout  this 
paper).  Eq.  (2)  is  derived  from  the  assumption  that  the  deviatoric  stress  is 
linearly  related  to  the  deviatoric  strain-rate  and  the  effective  deviator  is 
strain-rate  (D  is  related  to  the  effective  stress  iag)  as 


f  = 


n 
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2  £&j 


.  r  ,2  .  .  ,1/2 

where  r  =  (■=  e  .  .e  .  . ) 

3  ij  ij 


Noting  that  for  the  linear  viscous  material  the 
viscosity  is  obtained  as  a  slope  of  o-e  curve  and  s^..  eq.  (7)  becomes  equal 
to  2u*e^j,  we  can  set  jj  as 


H  = 


34 


(4) 


Here  ^s  assumed  to  be  related  to  the  volume  fraction  of  the  voids  f  in 
accordance  with  u  in  eq.  (1)  of  Appendix  A.  Then  the  total  stress  can  be 
written  by 


M 


1-n 


or.  . 


ij  =  ~Psij  ~  2m  etj 

where  p  is  the  hydrostatic  pressure  and  o^oQ  is  given  by 

c  1  7n 

(«! 

O'  rr_ 


(5) 


(6) 


The  extreme  cases  of  these  materials  are  a  linear  viscous  material  with  n  =  1 
and  a  rigid/plastic  material  with  n~>«. 


r, 
»'*  « 
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Let  us  introduce  a  representative  void  ft  of  radius  R  in  the  origin  of  the 
compressible  body  (Fig.  14)  and  apply  the  variational  principle  proposed  by 
Hill  C193  and  also  used  by  Glennie  C173  to  the  problem  of  dynamic  void  growth 
in  a  non-linear  compressible  viscous  material  D.  Following  Glennie,  the 
approximate  velocity  field  u^.  can  the  solution  of  the  problem  provided  that 


du . 


Du. 


3o* 


f  "V  Ot  -  -  ax,  ,MV  5  •  4aij"jnidS 


(7) 


for  any  velocity  field  w^(  where  and  are  the  stress  field  corresponding  to  tfie  ||| 

assumed  velocity  field  and  the  stress  field  at  the  outer  boundary,  *« 

V-  . 

respectively,  p  is  the  density  of  the  matrix  and  D/Dt  is  the  material  time  . 

derivative.  w 
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In  general,  the  velocity  can  be  approximated  as 


Qi  =  e°ijXj  ♦  Du?  +  Euf  (8) 

where  u^D  represents  a  spherically  symmetric  expansion,  while  denotes  the 
shape  distorsion  of  the  void.  Both  fields  should  vanish  at  large  distance  from 
the  void.  However,  noting  that  the  effect  of  the  shape  change  can  be  neglected 
under  high  stress  triaxiality  CIO, 203,  we  consider  only  the  dilational  growth 

field  u^D  an(j  U3e  the  3ame  velocity  field  as  those  used  by  Rice  and  Tracey  C203 
under  the  assumption  that  the  disturbance  velocity  field  is  limited  in  the 
small  vicinity  of  void  and  from  the  fact  that  the  pure  matrix  is 
incompressible.  Thus  we  have 


+  Du 


~D 


(9) 


where 
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=  (x^x^) 


1/2 


and  R  is  the  radius  of  a  representative  void.  The  normalizing  field  w^  in 
(7)  i3  defined  as 


wj  =  ^io) 

It  is  noted  here  that  u^0  an(j  w_.  die  out  a3  the  distance  r  from  the  origin 
increases.  After  some  manipulation  (see  the  details  in  C263),  we  arrive  at  the 
following  growth  equation 


i .  ^ 

R  3  2n  a^oo  ae»  2 


(1+2^) 
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(l+2^2) 

l  R2®a2>} 


+  (n  -  1 > ( n  t  0-432),  n 
n2 


(11) 


where  the  stress  ratio  (cj^/c^0)  for  various  boundary  conditions  must  be  given. 
The  effective  stress  at  far  field  ao°°  is  given  by 


ms  n  vttvw  vtt  rriivn  vivfw 


2M(l-?)eA 

£ - -} 

<*o 


1/n 


(12) 


However,  the  mean  stress  am'M  at  remote  boundary  cannot  be  given  explicitly 
except  for  the  linear  viscous  material  (n  =  1)  which  is  given  by 


°m  =K  ekk  '  £kk  =  (1  +  25)eA  <13) 

Hence,  the  stress  ratio  <x  for  n  =  1  is  set  to 

a  -  q»-  _  (-VL_t  m  _  (l^.v).  a  ±  m  (14) 

°>  2u  (1  -  E,)  3(1  -  2v)  (1  -  E,) 

and  it  is  assumed  that  eq.  (14)  holds  for  the  non-linear  material  in  order  to 
have  the  same  condition  on  a.  The  stress  ratio  a  in  eq.  (14)  is  dependent  not 
only  on  the  mode  parameter  E,  but  also  on  the  Poisson's  ratio  v,  thereby  a  is  a 
function  of  the  volume  fraction  of  voids  f. 

Next  let  us  consider  three  kinds  of  straining  mode,  i.e.,  ^  -  0,  -v  and  1. 
However,  when  the  initial  volume  fraction  of  voids  (fQ)  is  zero,  the 
surrounding  matrix  becomes  strictly  incompressible.  Then  the  stress  ratio  a 
given  by  eq.  (14)  becomes  infinite,  which  means  infinite  void  growth  except  for 
the  case  of  E,  =  -v.  When  E,  =  -v  (uniaxial  tension  mode),  «  becomes  finite 

(1/3).  Especially,  if  E,  -  -v  =  -1/2,  e^  =  o  as  seen  in  eq.  (13)  and  our 
growth  equation  can  give  rise  to  the  existing  solutions  in  which  the 
surrounding  matrix  is  assumed  to  be  incompressible.  Under  these  conditions  (E; 

=  -v  =  -1/2,  fQ  =  o),  the  effective  strain-rate  ep  is  equal  to  and  eq.  (11) 
is  reduced  to 
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By  setting  n  =  1  and  using  eq.  (3.6),  we  can  simplify  eq.  (3.28)  as 
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which  is  the  same  as  the  classical  solution  C183  derived  from  the  Bernoulli 
equation  except  that  the  last  terra  inside  the  bracket  is  added  and  pressure 
difference  is  replaced  by  here.  For  n  being  infinity,  we  have  from  eq. 


|  =  lira  -e-(RR+1R2-  —(e"  ,2)}+  ( )  ( n±Q-v.432)  * 
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which  coincides  with  the  Glennie's  solution  C173  by  using  o® 

e 

Furthermore,  if  we  neglect  the  inertia  term,  eq.  (15)  is  reduced  to 
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which  was  obtained  by  Budiansky  et  al  for  the  quasi-static  void  growth  in  a 
non-linear  viscous  material  C103.  Thus,  the  dynamic  void  growth  equation 
derived  here  can  include  the  existing  solutions  for  special  cases,  in  which  the 
material  is  assumed  to  be  incompressible. 

For  the  case  of  uniaxial  3train  mode  (£  =  0),  eq.  (11)  is  reduced  to 
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( n-1 ) ( n+0.432) 


where  «.  =  (l+v)/3(l-2v) .  Then  we  can  obtain  the  volume  fraction  of  voids  by 
solving  eq.  (19)  by  the  Runge-Kutta  method.  Since  we  will  simulate  the  void 
growth  by  the  present  model  in  a  ductile  metal  impacted  by  flyer  test  that 
gives  rise  to  the  state  of  uniaxial  strain  (Taya  et.  al  C223),  we  focus  only  on 
the  numerical  results  for  the  uniaxial  strain  mode  (£  =  0). 


Appendix  C 


The  Eshelby  s  tensors  are  given  in  C283.  We  present  them  in  a  more 
convenient  form  below.  The  suffices  ijk£  in  the  Eshelby 's  tensors  are 
consistent  with  the  coordiante  systems  (Fig. 14). 

1.  When  £1  is  a  sphere 
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2.  When  ft  is  prolate  spheroid  along  x3-axis 
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